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Abstract 

A theory is presented which allows us to accurately calculate the surface tensions and the 
surface potentials of electrolyte solutions. Both the ionic hydration and the polarizability are 
taken into account. We find a good correlation between the Jones-Dole viscosity fi-coefficient 
and the ionic hydration near the air-water interface. The kosmotropic anions such as fluoride, 
iodate, sulfate and carbonate, are found to be strongly hydrated and are repelled from the inter- 
face. The chaotropic anions such as perchlorate, iodide, chlorate and bromide are found to be 
significantly adsorbed to the interface. Chloride and bromate anions become weakly hydrated 
in the interfacial region. The sequence of surface tensions and surface potentials is found to 
follow the Hofmeister ordering. The theory, with only one adjustable parameter, quantitatively 
accounts for the surface tensions of 10 sodium salts for which there is experimental data. 
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Introduction 



Electrolyte solutions have been the subject of intense study for over a century. While the bulk 
properties of electrolytes are now quite well understood, their behavior at interfaces and surfaces 
still remains a puzzle. Over a hundred years ago, Hofmeister observed that presence of electrolyte 
in water modified significantly solubility of proteins. While some salts lead to protein precipitation 
(salting out), other salts stabilize proteins increasing their solubility (salting in). A few years after 
this curious observation, Heydweiller 1 discovered that salt dissolved in water increased the surface 
tension of solution-air interface. While cations had only a small influence on the surface tension, 
anions affected it quite significantly. Furthermore, the magnitude of the variation of the surface 
tension followed the same sequence discovered by Hofmeister earlier. It appeared that the two 
effects were related. 

Because of its great importance for biology, over the last century there has been a tremendous 
effort to understand the ionic specificity. Since air-water interface is relative simple, as compared 
to proteins, most of the theoretical work has concentrated on it. Langmuir 2 was the first to attempt 
a theoretical explanation of the physical mechanism behind the increase of the surface tension pro- 
duced by electrolytes. Using the Gibbs adsorption isotherm equation, Langmuir concluded that the 
phenomenon was a consequence of ion depletion near the air-water interface, and suggested that 
the depleted layer was about 4 A in width. No clear explanation for existence of this ion-free layer 
was provided by Langmuir, and soon it became clear that in order to obtain a reasonable agree- 
ment with experiments, its width had to be a function of ionic concentration. 3 A further insight 
was provided by Wagner, 4 based on the Debye-Hiickel (DH) theory of strong electrolytes. 5 Wag- 
ner argued that ionic depletion was a consequence of the electrostatic repulsion produced by the 
interaction of ions with their electrostatic images across the air-water interface. Wagner's theory 
was quite complicated and a simplified version was proposed by Onsager and Samaras (OS), who 
also derived a limiting law which, they argued, had a universal validity for all electrolytes at suf- 
ficiently low concentrations. 6 Indeed, careful experimental measurements have confirmed the OS 
limiting law. 7 However at larger concentrations, OS theory was found to strongly underestimate 
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surface tensions. 

Similar to Langmuir, Wagner and OS integrated the Gibbs adsorption isotherm equation to 
obtain the excess surface tension. Some time ago, Levin and Flores-Mena (LFM) proposed a dif- 
ferent approach based on the direct Helmholtz free energy calculation. 8 The LFM theory combined 
Langmuir, Wagner, and OS insights into one theory. They argued that in addition to the ion-image 
interaction, ionic hydration leads to a hard-core-like repulsion from the Gibbs dividing surface. 
Using a 2 A hydrated radius of Na + and CI", they were able to obtain a very good agreement 
with the experimental measurements of surface tension of NaCl solution, up to 1 M concentration. 
However, the LFM theory failed to correctly account for the ionic specificity, predicting that the 
surface tension of Nal salt should be larger than that of NaCl, contrary to experiments. It was 
clear that the LFM theory was still lacking an important ingredient. An indication of the missing 
ingredient was already present in the work of Frumkin, 80 years earlier. 9 Frumkin measured the 
electrostatic potential difference across the solution-air interface and found that for all halogen 
salts, except fluoride, potential was lower in air than in water. This meant that anions were pref- 
erentially solvated near the interface. Bostrom et al.— suggested that this ionic specificity, was a 
consequence of dispersion forces arising from finite frequency electromagnetic fluctuations. Their 
theory, however, predicted surface potentials of opposite sign to the ones measured by Frumkin, 
implying that cations were preferentially adsorbed to the interface. This was clearly contradicted 
by the simulations on small water clusters, 11 as well as by the subsequent large scale polarizable 
force fields simulations 12 and by the photoelectron emission experiments. 13 All these agreed with 
Frumkin that large anions, and not small cations, that are preferentially solvated near the air-water 
interface. A theoretical explanation for this behavior was advanced by Levin, 14 who argued that 
anionic surface solvation was a consequence of the competition between the cavitational and the 
electrostatic energies. The cavitational energy arises from the perturbation to the hydrogen bond 
network produced by ionic solvation. This results in a short range force which drives ions towards 
the interface. This is counterbalanced by the electrostatic Born solvation force which arises from 
the dipolar screening of the ionic self-energy in aqueous environment. For hard (weakly polariz- 
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able) ions, the Born energy is much larger than the cavitational energy, favoring the bulk solvation. 
However for large polarizable ions, the energy balance is shifted. For such ions, ionic charge can 
easily redistribute itself so that even if a large fraction of ionic volume is exposed to air, the elec- 
trostatic energy penalty for this remains small, since most of the ionic charge remains hydrated. 
This means that through surface solvation, large, strongly polarizable ions, can have the best of two 
worlds — gain the cavitational energy at a small price in electrostatic self energy. Levin derived 
the interaction potential quantifying this effect.— In a follow up paper, Levin et al.— used this 
potential to quantitatively account for the surface tensions of all sodium-halide salts. In this paper 
we will extend the theory of reference 15 to calculate the surface tensions and the surface potentials 
of others sodium salts and to derive the Hofmeister series. 



Consider an electrolyte solution confined to a mesoscopic drop of water of radius R, which cor- 
responds to the position of the Gibbs dividing surface (GDS).— ^ We define the adsorption (ion 
excess per unit area) as 



where p±(r) are the ionic density profiles and q, = p+(0) = P-(0) is the bulk concentration of 
electrolyte. If N ion pairs are inside the drop, eqQ] simplifies to T± = N /AtzR 1 — CbR/3. 

The water and air will be treated as uniform dielectrics of permittivities e w = 80 and £ a = 1, 
respectively. The surface tension can be obtained by integrating the Gibbs adsorption isotherm 
equation, dy = — r + d/i + — r dp , where /3/l± = ln(o,A±) are the chemical potentials and A± 
are the de Broglie thermal wavelengths. Let us first consider alkali-metal cations, such as lithium, 
sodium, or potassium. Because these cations are small, they have large surface charge density, 
which leads to strong interaction with surrounding water molecules, resulting in an effective hy- 
drated radius a/,. We can, therefore, model these ions as hard spheres of radius with a point 



Model and theory 




(i) 
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charge q located at the origin. Because of their strong hydration, these cations can not move across 
the GDS since this would require them to shed their solvation sheath. For mesoscopic drops we 
can neglect the curvature of the GDS. To bring a cation from bulk electrolyte to some distance 
z> from the GDS then requires-^' 



= l_ r dke -2s {z - ah) k[scoMkak) -k**h{ka h 
2e w Jo s[scosh(kah) +ks'mh(kah)\ 



(2) 



of work. The GDS is located at z = and the axis is oriented into the drop. We have defined s = 
\J K 2 + k 2 , where K = ^%nq 2 Cb/ £ w ksT is the inverse Debye length. The eq[2]is well approximated 
by 



W ap (z;a h ) 



W(a h ;a h )a h 



-2K{z—ah) 



(3) 



see Figure 1 This form will be used later to speed up the numerical calculations. 




z[A] 



Figure 1: Comparison between the complete expression (eqO and its approximation (eqO, for an 
ion of radius 2.5 Aat 1M concentration. 



Unlike small cation, large halogen anions of bare radius ciq have low electronic charge density 
and are weakly hydrated. The polarizability of an anion is y a and we define its relative polarizability 
as a = Ya/ciQ. We can model these ions as imperfect spherical conductors. When an anion is far 
from the interface, its charge — q is uniformly distributed over its surface. However, when an anion 
begins to cross the GDS, its charge starts to redistribute itself on the surface so as to leave most of 
it in the high dielectric environment. 14 The fraction of charge x which remains hydrated when the 



ionic center is at distance z from the GDS is determined by the minimization of the polarization 
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energy, 
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where 6(z) = arccos[— z/ao] and g = (1 — a) /a. We find 
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Substituting this back into eq |4l we obtain the polarization potential. This potential is repulsive, 
favoring ions to move towards the bulk. Nevertheless, the repulsion is quite soft compared to 
a hard-core-like repulsion of strongly hydrated cations. The force that drives anions toward the 
interface arises from cavitation. When ion is dissolved in water it creates a cavity from which 
water molecules are expelled. This leads to a perturbation to the hydrogen bond network and an 
energetic cost. For small voids, the cavitational energy scales with their volume. 17 As an ion moves 
across the interface, the cavity that it creates in water diminishes proportionally to the fraction of 
the volume exposed to air,— ^ producing a short-range interaction potential that forces ions to 
move across the GDS, 



U r . 



vciq for z > ciq , 



>o(i + 1 ) ( 2 ~i) ^-a <z< 



(6) 



a 



where v ~ 0.3ksT / A 3 is obtained from bulk simulations. 18 For small strongly hydrated cations, 
this energy gain does not compensate for the electrostatic energy penalty of exposing ionic charge 
to the low dielectric environment. For large polarizable ions, on the other hand, electrostatic energy 
penalty is small, and the cavitational energy is sufficient to favor the surface solvation. The total 
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potential felt by an unhydrated anion is then 
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U tot (z) = { 



W(z;ao) + val+2^ for z>a , 
W(a ;a )z/ao + U p (z) + U cav (z) for < z < a 
U p (z) + U cav (z) for -a <z<0. 



(7) 



The density profiles can now be calculated by integrating the non-linear modified Poisson- 
Boltzmann equation (mPB): 



V 2 0(r) 

p+W 

p_(r) 



Aizq 



[p+W-p-W] 



Ne M(r)-pU tot (r) 



(8) 



where © is the Heaviside step function. To speed up the numerical calculations we can replace 



ap- 



Sodium-Halogen Salts 

First, we study the sodium-halogen salts. 15 The anion radii were obtained by Latimer, Pitzer and 
Slansky 19 by fitting the experimentally measured free energies of hydration to the Born model. 
Since our theory in the bulk also reduces to the Born model, these radii are particularly appropriate: 
ai = 2.26 A, dBr = 2.05 A, aci = 1.91 A and clf = 1.46 A. For ionic polarizabilities we use the 
values from the reference:- 2 ^ y 7 = 7.4 A 3 , y Br = 5.07 A 3 , y a = 3.77 A 3 and y F = 1.31 A 3 . 

The excess surface tension can be obtained by integrating the Gibbs adsorption isotherm equa- 
tion (eqQ])- We start with Nal. Since I~ is large and soft, it should be unhydrated in the interfa- 
cial region. Adjusting the hydrated radius of Na + to best fit the experimental data— for Nal, we 
obtain 15 ah = 2.5 A. We will use this partially hydrated radius of Na + in the rest of the paper. 
Considering that Br" is also large and soft, we expect that it will also remain unhydrated in the 
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Figure 2: Surface tensions for NaF, NaCl, NaBr and Nal. The symbols are the experimental 
data 21 " 23 and the lines are the results of the present theory. The inset shows the surface potential 
difference as a function of molar concentration. 



interfacial region. This expectation is well justified and we obtain a very good agreement with the 
experimental data, 23 see|Figure~2| For F" the situation should be very different. This ion is small, 



hard and strongly hydrated. This means that just like for a cation, a hard core repulsion from the 
GDS must be explicitly included in the mPB equation. For hydrated (or partially hydrated) anions 
the density is then 



N®(R -a h - r ) e M(r)-pW(z;a h ) 
$*~ ah Anr 2 dreP^ r )-P w & a h) 



(9) 



Using this in eq [81 an almost perfect agreement with the experimental data 21 is found for NaF, 
using the usual bulk hydrated radius of F~, = 3.52 A,— see |Figure 2| The table salt, NaCl, 



is the most difficult case to study theoretically, since Cl~ is sufficiently small to remain partially 
hydrated near the GDS. We find a good fit to the experimental data, 22 |Figure 2[ using a partially 
hydrated radius of Cl~, ah = 2.0 A, which is very close to its bare size. 



The ionic density profiles for these salts are plotted in Figure 3 In agreement with the po- 
larizable force fields simulations, 11112 the large halogens I~ and Br~ are strongly adsorbed at the 
interface. However, just as was found using the vibrational sum-frequency spectroscopy, 25 their 
concentration at the GDS remains below that in the bulk. The electrostatic potential difference 
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Figure 3: Theoretical ionic density profiles for NaF, NaCl, NaBr and Nal at 1M. The GDS is at 
r = 300 A. 



across the air-water interface is % = — 0(0). We find that our theoretical values for the sur- 
face potentials are in reasonable agreement with the measurements of Frumkin 9 - 26 and Jarvis et 
al.-^ ITable l l 

Table 1: Surface potentials difference at 1 M for various salts 

Frumkin 9 - 26 (mV) 



NaF 



NaCl 



NaBr 



Nal 



NalO, 



NaBrO, 



NaNO, 



NaCIO, 



NaCIO, 



Na 2 CQ 3 
Na 2 SQ 4 



Calculated (mV) 



4.7 



-2.1 



-9.4 



■14.3 



-0.12 



-8.27 



■11.02 



-31.1 



10.54 
10.17 



■1 



-39 



-17 



-41 



-57 



Jarvis et al. 27 (mV) 



-1 



-21 



35 



Oxy-Anion Salts 

We next study the surface potentials and the surface tensions of sodium salts with more complex 
anions, such as: chlorate, nitrate, bromate, iodate, perchlorate, sulfate and carbonate. Naively 
one might expect that because of their large size, these oxy-anions will loose their hydration near 



the GDS. This, however, is not necessarily the case. At the moment there is no reliable theory 
of hydration, nevertheless, it has been known for a long time that ions can be divided into two 
categories: structure-makers (kosmotropes) and structure-breakers (chaotropes). This dichotomy 
is reflected, for example, in the Jones-Dole viscosity 5-coefficient. It is found experimentally that 
for large dilutions, viscosity varies with the concentration of electrolyte as:~ 

7] r = l+Ap l / 2 +Bp , (10) 

where A is a positive constant, resulting from the ion-ion interactions; and B is related to the ion- 
solvent interactions. For kosmotropes .^-coefficient is positive, while for chaotropes it is negative, 
see ITable~2l In particular, we observe that F is a kosmotrope, while Br" and I - are chaotropes. 
Chloride ion appears to be on the border line between the two regimes. The classification of ions 
into kosmotropes and chaotropes correlates well with our theory of surface tensions of halogen 
salts — kosmotropes were found to remain strongly hydrated near the GDS, while chaotropes 
lost their hydration sheath near the interface. It is curious that the bulk dynamics of ion-water 
interaction, measured by the value of the viscosity 5-coefficient, is so strongly correlated with the 
statics of ionic hydration, measured by the surface tension of electrolyte solutions. 



Table 2: 5-coefficients for various ions.— 



Ions 


5-coefficient 


Ions 


5-coefficient 


Na + 


0.085 


Br0 3 


0.009 


F" 


0.107 


N0 3 


-0.043 


cr 


-0.005 


CIO3 


-0.022 


Br" 


-0.033 


C10 4 


-0.058 


r 


-0.073 


CO3 2 - 


0.294 


io 3 " 


0.140 


S0 4 2 " 


0.206 



It is reasonable, then, to suppose that the classification of ions into chaotropes and kosmotropes, 
based on their 5-coefficient, will also extend to more complicated ions as well. Thus, we expect 
that chaotropes NO3 , CIO3 and CIO4, will loose their hydration sheath near the GDS, while 
the kosmotropes IO3 , COf and SO^ - , will remain strongly hydrated. Furthermore, since B- 
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coefficients of these ions are large (as compared to F - ), we expect that these ions will remain as 
hydrated near the GDS as they were in the bulk, similar to what was found for fluoride anion. On 
the other hand, the .B-coefficient of BrOj is very close to zero and, similarly to CI - , we expect 
bromate to be only very weakly hydrated near the GDS. 

There is an additional difficulty with studying oxy-anions. Since these ions are not spherically 
symmetric, their radius is not well defined. Nevertheless, it has been observed that empirical radius 

00 = ^(^ + 1.4 A) (11) 

where d is the M-0 covalent bond length in the corresponding salt crystal and n oxy is the number 
of oxygens in anion, correlates very well with the experimental measurements of entropies of 
hydration. 29 Using this formula we calculate: a^o 3 = 1-98 A, ac/o 3 = 2.16 A and cicio 4 = 2.83 A, 
for the bare radius of the chaotropic oxy-anions. The polarizabilities of NO3 and CIO4 are given in 
the reference: 20 Jno^ = 4.09 A 3 and JqiOa = 5.4 A 3 . Unfortunately, this reference does not provide 
the polarizability of chlorate ion. However, since the polarizability is, in general, proportional to 




concentration [M] 

Figure 4: Surface tensions for NaI0 3 , NaBr0 3 and NaC10 3 . The symbols are experimental data 23 
and the lines are obtained using the present theory. The inset shows the electrostatic surface poten- 
tial difference and has the same horizontal axis label as the main figure. 

the volume of an ion, we can easily estimate it based on the polarizability of iodate, which is given 
in the reference, 20 Yio 3 = 8-0 A 3 . We then obtain Ycio 3 = 5-3 A 3 . Using eq[8l we can now calculate 
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the surface potential, the ionic density distribution and, integrating the Gibbs adsorption isotherm, 



the excess surface tension, Figure 4, Figure 5, Figure 6 t Figure 7, and Figure 8 We stress that in 
all the calculations we are using the same value of sodium hydrated radius, = 2.5 A, obtained 
for halogen salts. For nitrate and chlorate, the predictions of the theory are in good agreement with 



the experimental measurements of Matubayasi, 23 see |Figure~4] and [Figure 5 




0.2 
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Figure 5: Surface tension of NaN0 3 . The symbols are the experimental data 23 and the lines are 
the present theory. The inset shows the electrostatic surface potential difference vs. the molar 
concentration. 
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Figure 6: Surface tensions of NaC10 4 . The symbols are the experimental data 23 and the lines are 
the present theory. The solid line is obtained using acw 4 = 2.83 A from eq. (PTTI) . and the dashed 
line using acio 4 — 2.75 A. The inset shows the electrostatic surface potential difference vs. the 
molar concentration. 



For perchlorate, the theoretical surface tension is slightly lower than what was found experi- 



12 



mentally, |Figure 6| The difficulty is that since the cavitational energy scales with the volume of the 



ion, for large ions it becomes very sensitive to the precise value of the ionic radius. For example, 
if we use a slightly smaller radius of perchlorate acio 4 = 2.75 A, we obtain a good agreement with 



experiment, Figure 6 The theory shows that C10 4 is adsorbed stronger than either C10 3 or I , 
Figure 7|and|Figure 8| This is conclusion is also in agreement with the simulations.— 



see 
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Figure 7: Theoretical ionic density profiles for NaI0 3 , NaBr0 3 , NaN0 3 and NaC10 3 at 1M. The 
GDS is at r = 300 A. 
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Figure 8: Theoretical ionic density profiles for Na 2 C0 3 , Na 2 S0 4 and NaC10 4 at 1M. The GDS is 
at r = 300 A. 



We next calculate the surface tension of sodium salts with kosmotropic anions. Since the B- 
coefficient of iodate is so large, we expect that this ion will remain completely hydrated near the 
GDS. The bulk hydrated radius of I0 3 is ah = 3.74 A. 24 Using this in eq|H we obtain an excellent 
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agreement with the experimental data, 23 |Figure 4[ The 5-coefficient of bromate is almost zero, 



so that Br0 3 should be only very weakly hydrated near the GDS. Indeed, using a/, = 2.41 A, we 
obtain a very good agreement with the experimental data, 23 |Figure 4| This value is only slightly 



above the crystallographic radius of Br0 3 , agrO^ = 2.31 A, obtained using eqfTTl 
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Figure 9: Surface tensions of Na 2 C0 3 and Na 2 S0 4 . The circles and the squares are the experimen- 
tal data for sodium carbonate and soldium sulfate, respectively 21 . The up and down triangles are 
the data of Jarvis and Scheiman. 27 The plus symbols are the data of Weissenborn and Pugh 31 for 
the Na 2 S0 4 salt. The lines are the present theory. The inset shows the surface potential difference 
vs. the molar concentration. 



Finally, we consider the surface tension of sodium salts with divalent anion. The predictions 
of the theory, in this case, should be taken with a grain of salt, since it is well known that Poisson- 
Boltzmann theory begins to loose its validity for divalent ion as a result of inter-ionic correla- 
tions. 32 Nevertheless, it is curious to see how well the theory does in this limiting case. The 

9 

5-coefficients of SO4 and C0 3 are large and positive, so that we expect these ions to be fully 
hydrated, a/, = 3.79 A and 3.94 A, 24 respectively. Calculating the surface tensions, we find that our 
results are somewhat above those measured by Matubayasi et al.— . Curiously, the theory agrees 
well the measurements of Jarvis and Scheiman 27 and of Weissenborn and Pugh, 31 Figure 9 At 
this point we do not know what to make of this agreement. The surface potentials of sodium sulfate 
and sodium carbonate are listed in lTable fl In ITable 31 we summarize the classification of all ions 
and the radii used in calculations. 
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Table 3: Ion classification into chaotropes (c) and kosmotropes (k) and their effective radii. The 
kosmotropic radii listed are either hydrated or partially hydrated. The chaotropic radii are the bare 
ones, for oxy-anions given by eqfTTl 



Ions 


chao/kosmo 


; 

radius ( A") 


Na + 


k 


2.5 


F" 


k 


3.54 


cr 


k 


2 


Br~ 


c 


2.05 


I- 


c 


2.26 


IO3 


k 


3.74 


Br0 3 


k 


2.41 


N0 3 


c 


1.98 


CIO3 


c 


2.16 


CIO4 


c 


2.83 


C0 3 2 


k 


3.94 


S0 4 2 " 


k 


3.79 



Conclusions 

We have presented a theory which allows us to quantitatively calculate the surface tensions and 
the surface potentials of 10 sodium salts with only one adjustable parameter, the partial hydrated 
radius of the sodium cation. We find that anions near the GDS can be classified into kosmotropes 
and chaotropes. Kosmotropes remain hydrated near the interface, while the chaotropes loose their 
hydration sheath. For all salts studied in this paper, the classification of ions into the structure mak- 
ers/breakers agrees with the bulk characterization based on the Jones-Dole viscosity 5-coefficient. 
If anions are arranged in the order of increasing electrostatic surface potential difference, we obtain 
the extended lyotropic (Hofmeister) series, C0 3 2 - > S0 4 2 " > I0 3 > F" > Br0 3 > Cr > NO3 > Br~ 
> CIO3 > I" > CIO4 . To our knowledge this is the first time that this sequence has been derived 
theoretically. The theory also helps to understand why Hofmeister series is relevant to biology.— 
In general, proteins have both hydrophobic and hydrophilic moieties. If the hydrophobic region 
is sufficiently large, it can become completely dewetted, 17 surrounded by a vapor-like film into 
which chaotropic anions can become adsorbed. This will help to solvate proteins, but at the same 
time will destroy their tertiary structure, denaturing them in the process. On the other hand, the 
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kosmotropic ions will favor protein precipitation. The action is again two fold. The kosmotropes 
adsorb water which could, otherwise, be used to solvate hydrophilic moieties. They also screen the 
electrostatic repulsion, allowing charged proteins to come into a close contact and to stick together 
as the result of their van der Waals and hydrophobic interactions. This leads to formation of large 
clusters which can no longer be dispersed in water and must precipitate. 
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